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NOTATIONS  AND  NUMERICAL  VALUE  USED 


a 

b 

r'6 

a 

<P 

A,B,C,D 

u 

d 

G 


v 

T 


Ta>Tt 


E 


a 

P 


inner  radius,  1" 
outer  radius,  2" 
cylindrical  coordinates 
normal  stress 
yield  stress,  170  ksi 
Airy  stress  function 
superposition  constants 
displacement 

coefficient  of  dislocation 

shear  modulus 

Poisson's  ratio,  0.3 

thermoelastic  potential 

temperature  at  r 

temperature  at  r=a,  r=b 

Young's  modulus,  30x10*  psi 

coefficient  of  thermal  expansion,  6.8x10 

radius  of  the  autofrettaged  interface 

temperature  at  r=p 
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FULLY  AUTOFRETTAGED  CASE 

The  plane  strain  stress  distribution  of  a  fully  autofrettaged 
tube  using  von  Mises  yield  condition  and  the  incompressibility  condi¬ 
tion  is  given  by 

20n  „  v  n  2 

d-~)logi}  (1) 


0r  =  ~  {log  -  -  —~~ 
v^3  b  b2-a2 


h2M  b  i 

_)  log  -  } 
rz  a 


ae  • 


20Q  y  _  2 

—  (1  ♦  log  L  -  a 


b2-az 


(1  ♦  —jO log  -  } 
rz  a 


/3  b 

This  distribution  can  be  simulated  either  by  a  dislocation.  Figure  1, 
or  by  a  steady  state  thermal  loading. 


(2) 


Figure  1.  A  portion  of  the  ring  between  two  adjacent  cross  sections 
is  cut  out.  If  the  ends  of  the  ring  are  joined  again, 
stresses  thus  produced  may  simulate  the  residual  stresses 
due  to  autofrettage. 
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Dislocation  Solution; 

Using  biharmonic  Airy  stress  function*  the  dislocation  solution 
can  be  obtained  by 

4>  ■  A  log  r  +  Br2  +  Cr2  log  r  (3) 

The  dislocation  is  expressed  by  the  jump  condition 

0=2tt 

[2Gu0]  «=  d  •  r  (4) 

0=0 


This  condition  together  with  traction  free  conditions  at  the  inner  and 
outer  radii  gives 


_i_  i!>>L  :„*  >> 

4ir(l-\>)  b2-a2  a 


- ^ -  {-I-* log  —  +  1  +  2  log  b} 

16ir(l-v)  (b2-a2)  a 


Using  the  formulas 


C 


d 

8ir(l-v) 


1  324»  ,  1  34> 
r2  302  r  3r 

3ad> 

3r2 


(5) 


(6) 


T 


Timoshenko,  S.  and  Goodier,  J.  N. ,  Theory  of  Elasticity.  McGraw-Hill 
Co.,  1951,  2nd  Edition,  p.  56. 


The  stress  distribution  is  then  obtained  from  (3)  and  (5)  as 


0  =  - ^ - {log  L  - 

4tt(1-v)  b  b2-a2 

d 


a2  (1  -  £l)  log  4} 
r2  a 


00  = 


>2  b2 


{log  £  +  i  -  -4 — r  (i  +  4_0  log  4} 

4ir(l-v)  b  b  -az  rz  a 


The  equivalence  between  (7),  (8)  and  {1),  (2)  is  easily  seen  with 
dislocation  and  yield  stress  related  by 

H  20/s 


4tt(1-v)  /7” 


(7) 

(8) 

(9) 


Solution  of  Thermal  Loading: 

Using  the  superposition  of  Airy  stress  function  4>  and  thermo¬ 
elastic  potential  ip  (ref.  2),  the  solution  can  be  symbolically  written 
as 


[S]  =  A:[>-r2]  +  B^-rhogr]  +  Cjf^-logr]  +  D^cfr-r2]  (10) 
with  Ta  and  Tjj  as  steady  state  temperatures  at  the  inner  and  outer 
radii  respectively  and  using  the  traction  free  boundary  conditions  we 
have 

a  _  Ea  r-r  .  (1+loga)  (Tg-Tb)  1 

A1  '  1^7  LTa  log(b/a)  J 

B  -  "  Ea  (Ta-Tt>) 

1  4 ( 1 — v)  log (b/ a) 

,a2.2  (ID 

C1  =  -  B1  l0^b/a) 

D1  =  A1  +  1  B1  [1  +  21ogb  +  4^3-  log  (4)] 

A  *  b  -a  a 


2Sadowsky,  M.  A.  and  Hussain,  M.  A.,  "Thermal  Stress  Discontinuities 
in  Microfibers, "  Watervliet  Arsenal  Technical  Report  WVT-RR-6401, 
April  1964. 
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Using  the  formulas 


1  lit .  IM 

r5"  36*  r  3r 


Je 


_  lit 


(12) 


3r2 


The  stress  distribution  is  obtained  from  (6),  (10),  (11)  as 

b 


°r  = 


'0 


_ Ect(Ta.:Ti) _ dog  i  -  sL 

2(l-v) log(b/a)  b  b2-a2 


(l  -  li)log  °) 
r^  a 


(13) 


EaDV-Tb) 


(1  +  log  - 


(l+li)log£}  (14) 
r2  a 


2(l-v) log (b/a)  b  b2-a2 

The  equivalence  between  (13),  (14)  and  (1),  (2)  is  easily  seen  with 
the  temperature  gradient  and  yield  stress  related  by 

E<*(Ta-Tb)  2oq 


(15) 


2(l-v)log  (b/a)  /T 
PARTIALLY  AUTOFRETTAGED  CASE 

The  plane  strain  stress  distribution  of  a  partially  autofrettaged 
tube,  using  the  same  von  Mises  and  incompressibility  conditions  as 
before,  is 


H  {(21°*  §  •  1  ♦ 

-  Pl<pr- 

a<r<p 

(16) 

p£r<b 

(17) 

<(2lo»  f  ♦ 1  * 

-  p  (_L  +  JL)> 

1  b2  r2 n 

a£r<p 

(18) 

p<r<b 

(19) 

4 


where  ?i  =  -a0f  [d  -  ♦  21og(p/a)]. 


Solutions  of  Thermal  Loading: 

Using  the  superposition  of  Airy  stress  function  <f>  ar.d  thermo- 
elastic  potential  \p,  it  is  sufficient  to  write  the  solution  symbolically 
as 


[S] 


A2[^-r2]  +  B2[^-r2logr]  +  C2[^-logr]  +  D2 [cp-r2 ]  ,  a£r<p  (20) 


In  order 


>A3[>-r2]  +  C3[ip-logr]  +  D3[p-r2]  ,  p<r<b  (21) 

to  obtain  stress  distribution  given  by  (16) -(19)  we  must  have 


^2  -  °2  4  [2  +  21°gP  “  uT  fP2-pi)l 


A 


B,  =  - 


7T 


C2  =  "  pi  A 


(22) 


a3  -  d3  -  -  a 

°3  =  £ 


The  temperature  profile  from  (20)  and  (21)  is 

4A-  +  4B~(1  +  logr)  ,  a<r<p 

(23) 

4A3  p_<r<b 

It  is  seen  that  the  temperature  is  constant  in  the  outer  region, 
p<i<b,  and  logarithmically  distributed  in  the  inner  region,  a.<r<p. 

Let  Ta,  Tp  be  the  temperatures  at  r  =  a,  and  r  =  p  respectively.  These 
temperature  boundary  conditions  give  the  equivalence  between  the 
temperature  gradient  and  the  yield  stress 


EctT  = 
(1-v) 
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Eot  (Ta-Tp)  2a0 
2(l-v)log(p/a)  /j 


(24) 


The  temperature  profile  of  (23)  is  then  given  by 


T  =  T  -  iTrM  log(r/a)  a<r<p 
a  log (p/a) 


(25) 


T  =  Tp  p<r<b 

Once  the  temperature  distribution  is  known,  all  the  remaining  super¬ 
position  constants  can  be  specifically  determined.  It  should  be  noted 
that  we  have  neglected  the  axial  stress  computation  which  can  easily 
be  taken  care  of  by  the  method  discussed  on  page  409  of  Reference  1. 


A  NUMERICAL  EXAMPLE 

Consider  a  tube  of  inner  radius  a  *  1,  outer  radius  b  *  2,  with 
material  constants  E  -  30x10*  psi,  v  *  0,3,  a  *  6. 8x10" 6  in/in/°F, 
a0  ■  170x10*  psi;  the  temperature  distribution  was  computed  from  (25) 
for  30%,  60%  and  100%  autofrettaged  cases,  shown  in  Figure  2.  Using 
these  temperature  distributions  as  temperature  input  in  a  finite 
difference  computer  program  based  on  the  theory  of  thermal  stress  in 
section  9-10  of  Reference  3  we  obtain  the  stress  distributions.  The 
results  are  compared  in  Table  I  with  the  exact  solution  given  by  (14), 
and  are  also  graphically  shown  in  Figure  3. 


^Timoshenko,  S.  and  Goodier,  J.  N.,  Theory  of  Elasticity,  McGraw-Hill 
Co.,  1951,  2nd  Edition,  p.  56. 

3Boley,  B.  A.  and  Weiner,  J.  H. ,  "Theory  of  Thermal  Stresses,"  John 
Wiley  6  Sons,  1960. 
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Percent  Overstrain 


Exact  Solution 


Finite  Difference 


30% 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 
1.9 
2.0 

-  92190 

-  48446 

-  12325 

18205 

16442 

15020 

13856 

12891 

12083 

11398 

10814 

-  92897 

-  49567 

-  13636 

16825 

15236 

13948 

12890 

12010 

11271 

10643 

10106 

1.0 

-143955 

-145316 

1.1 

-  95719 

-  97552 

1.2 

-  56182 

-  58225 

1.3 

-  22993 

-  25102 

1.4 

5422 

3323 

60% 

1.5 

30153 

28107 

1.6 

51978 

50007 

1.7 

48359 

46593 

1.8 

45326 

43724 

1.9 

42759 

41289 

2.0 

40568 

39205 

1.0 

-166539 

-168854 

1.1 

-116343 

-119099 

1.2 

-  75316 

-  78246 

1.3 

-  40966 

-  43929 

1.4 

-  11631 

-  14550 

100% 

1.5 

13842 

11005 

1.6 

36275 

33539 

1.7 

56267 

53640 

1.8 

74269 

71749 

1.9 

90621 

88206 

2.0 

105590 

103274 

Figure  3.  Thermal  stresses  obtained  from  Eq.  (14)  using  temperature 
distributions  shown  in  Figure  2  simulating  30%,  60%,  and 
100%  overstrain. 
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The  same  temperature  distribution  was  also  used  in  the  finite 
element  NASTRAN  program  using  a  TEMP  (LOAD)  request  in  the  case  control 
deck.  The  results  are  again  compared  with  the  exact  solution  in 
Table  II. 

CONCLUSION 

A  simple  method  has  been  devised  to  simulate  partial  autofrettage 
residual  stresses  in  thick  walled  cylinders. 


TABLE  II.  COMPARISON  OF  aQ  (PSI)  WITH  FINITE  ELEMENTS 


Percent  Overstrain 

r 

Exact  Solution 

Finite  Element  fNASTRAN) 

m 

-  80392 

-  80638 

-  38795 

-  38981 

• 

■m 

-  4231 

-  4375 

1.325 

17727 

17818 

1.425 

16058 

16143 

30% 

1.525 

14707 

14782 

1.625 

13597 

13671 

1.725 

12675 

12740 

1.825 

11900 

11964 

1.925 

11243 

11307 

mm 

-130910 

-131102 

-  85128 

-  85261 

-  47361 

-  47454 

1.325 

-  15490 

-  15558 

1.425 

11915 

11872 

60% 

1.525 

35856 

35827 

1.625 

51010 

51136 

1.72S 

47551 

47660 

1.825 

44645 

44754 

1.925 

42179 

42283 

wm 

-152950 

-153095 

■ 

-105342 

-105430 

■ 

-  66178 

-  66228 

1.325 

-  33215 

-  33240 

1.425 

-  4938 

-  4941 

100% 

1.525 

19709 

19716 

1.625 

41482 

41496 

1.725 

60939 

60972 

1.825 

78500 

78538 

1.925 

94484 

L  - 

94519 

1.  Timoshenko,  S. 
Hill  Co.,  1951 

2.  '  Sadowsky,  M,  A, 

in  Microfibers 
April  1964. 
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DRDAR-LC B-TL,  WATERVLIET  ARSENAL,  WATERVLIET,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST  (CONT) 


NO.  OF 
COPIES 


NO.  OF 
COPIES 


COMMANDER 

DEFENSE  TECHNICAL  INFO  CENTER 
ATTN:  DTIA-TCA 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 

METAIS  &  CERAMICS  INFO  CEN 
BATTELLE  COLUMBUS  LAB 
505  KING  AVE 
COLUMBUS,  *10  43201 

MECHANICAL  PROPERTIES  DATA  CTR 
BATTELLE  COLUMBUS  LAB 
505  KING  AVE 
COLUMBUS,  *10  43201 

MATERIEL  SYSTEMS  ANALYSIS  ACTV 
ATTN:  DRXSY-MP 
ABERDEEN  PROVING  GROUND 
MARYLAND  21005 


COMMANDER 

US  ARMY  RESEARCH  OFFICE 
P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK,  NC  27709 
CCMMANDER 

US  ARMY  HARRY  DIAMOND  IAB 
ATTN:  TECH  LIB 
2800  POWDER  MILL  ROAD 
ADELPHIA,  MD  20783 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENG  ACT 

ATTN:  DRXPE-MT 

ROCK  ISIAND,  IL  61201 

CHIEF,  MATERIA IS  BRANCH 
US  ARMY  R&S  GROUP,  EUR 
BOX  65,  FPO  N.Y.  09510 

CCMMANDER 

NAVAL  SURFACE  WEAPONS  CEN 
ATTN:  CHIEF,  MAT  SCIENCE  DIV 
DAHLGREN,  VA  22448 

DIRECTOR 

US  NAVAL  RESEARCH  IAB 
ATTN:  DIR,  MECH  DIV 

CODE  26-27  (DOC  LIB) 
WASHINGTON,  D.  C.  20375 

NASA  SCIENTIFIC  &  TECH  INFO  FAC 
P.  0.  BOX  8757,  ATTN:  ACQ,  BR 

baltimore/washingtcn  intl  airport 

MARYIAND  21240 


NOTE:  PLEASE  NOTIFY  CCMflANDER,  ARRADCCM,  ATTN:  BENET  WEAPONS  LABORATORY, 
DRDAR-LCB-TL,  WATERVLIET  ARSENAL,  WATERVLIET,  N.Y.  12189,  OF  ANY 
REQUIRED  CHANGES. 


